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1 Introduction 



Many problems in particle physics and condensed matter physics can be mapped to anal- 
ogous problems in integrable quantum field theories and/or solvable lattice models in 
two dimensions by employing bold simplification or approximation. Understanding of 
integrable models has helped us gain insight into non-perturbative properties of four- 
dimensional quantum field theories, such as chiral symmetry breaking and confinement. 

In some physical systems, effects of boundaries are very important and hence for physi- 
cal application we have to extend integrable models to cases with boundaries. Methods for 
constructing integrable models with boundaries have recently been developed |^, ^, Q . 
One systematic approach lies in finding the refiection matrix (or K-matrix) obeying the 
boundary Yang-Baxter equation (BYBE) (or refiection equation). Only few models with 
integrable boundary terms in full generality have so far been constructed: (super) sine- 
Gordon theory |5|, massive Thirring model [^, affine Toda field theories 0, XX Z (and 
XXX) [|, I and XYZ spin-i chains flO . 



Quantum spin chains have played important roles in many aspects of physics such as 
strongly correlated electron systems. Recently there has been experimental interest in 
quasi one-dimensional systems of molecules which mimic spin chains of higher spin s. To 
our knowledge, however, for spin chain with s > 1 only boundary terms corresponding to 
diagonal K-matrices are known (see, e.g., [pT| ). 

In this paper we consider the 19- vertex model which is associated with the spin-1 
representation of Uq{sl2)- This model is the simplest extension of the six- vertex model. 
We construct and classify all refiection K-matrices of this model. From these K-matrices 
we derive the general integrable boundary terms for the XXZ spin-1 chain. As a special 
case we reproduce the hamiltonian obtained in |11|] by using the diagonal K-matrix. 



2 Nineteen-vertex Model 

The hamiltonian of the XXZ spin-1 chain is derived from the transfer matrix of the 
19-vertex model. This vertex model is defined in terms of the Boltzmann weights given 
by the R-matrix of spin-1 representation of Uq{sl2)- It is the trigonometric solution of the 
Yang-Baxter equation (YBE): 

(l ® R{u - u')) (i?(M) ® l) (l ® R{u')) = (/?(«') ® l) (l ® R{u)) (i?(u - u') O l) , (2.1) 
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where R{u) = PR{u) and P is the transposition matrix: P vi ® V2 = V2 ® vi. R{u) is a 
9x9 matrix with 19 non-zero entries [Q: 



R{u) = p{u) 



ai 



\ 



02 03 

0.4 05 dg 

03 O2 

O5 O7 05 

O2 03 
Og O5 04 

oa 02 



V 



01 = sin(M + 2?]) 

02 = sinu 
sin 27] 

sinw sin('U — rf) 



03 
04 



Oi ) 



05 



07 



sin(-u + Tj) 
sin 2r7 sin u 

sin(M + r]) 
sin 7] sin 27] 

sin(M + rj) 

og + sin-u, 



where p{u) ^ is an arbitrary function. If we choose p{u) = p±{u) 



(2.2) 



sin(2r7 ± u) ' 

then R satisfies 

R{u)R{-u) = 1, R{0) = 1. 

In the case of periodic boundary condition, the YBE implies a commuting family of 
transfer matrices and hence the integrability of the model. 



(2.3) 
(2.4) 



3 Reflection K-matrices 



We now consider the 19- vertex model with boundaries. As shown by Sklyanin [0], inte- 
grable models with boundaries can be constructed out of a pair of reflection K-matrices 
Kj-{u) in addition to R{u). K^{u) and K^{u) describe the effects of the presence of 
boundaries at the left and the right, respectively, and they both obey the BYBE. 
Let us construct the K-matrix which satisfies the BYBE: 



R{u-u')[k{u)®i)R{u+u)(i®K{u')) = {1(»K{u'))r{u+u)(^K{u)®i)R{u-u') 
In the present case, R{u) is given in (^.2|) and K{u) can be parametrized as 



/ xi{u) yi{u) z{u) \ 



K{u) 



p^{u) 



X2[U) 



y2[u) 



(3.1) 



(3.2) 



where p (u) ^ is an arbitrary function. 

The BYBE consists of 81 equations, many of which are dependent. To solve them, we 
proceed in a pragmatic fashion. The whole procedure is quite involved and requires a bit 
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of computer work. Here we only present the final result. More details of the derivation 
are relegated to the Appendix. 

We pick up a few easy-looking and independent equations. They are functional equa- 
tions and can be reduced to first order differential equations or algebraic equations. The 
solutions to these equations contain several arbitrary parameters. We then substitute 
these partial solutions back into the BYBE, which gives a set of equations with the num- 
ber of equations now less than 81 (since some equations are already satisfied). From this 
set of equations we pick up some simple-looking ones, solve them and determine the rela- 
tions among the parameters. Repeating this process for a couple of times, we end up with 
the solutions with three arbitrary parameters satisfying all 81 functional equations. The 
solutions obtained in this way exhaust all solutions of (|3.1| ) (for the 19- vertex model). 

All solutions can be classified into the following three cases: 

Class (A): sin2r/ ^ 0. 

yi {u) = /i sin(^ — 2 + -u) sin 2u, (m) = /i sin(^ — 2 + u) sin 2u, 

y2{u) = fism{( + ^ — u) sin2u, y2{u) = /isin(C + ^ — u) sin 2m, 

, , 2sin(?7-2M) ^ ^2^in{r]-2u) . ^ 
ziu) = fi sm2u, z[u) = jj sm2M, 

2 cos 7] 2 cos 7] 



Xi(u) = sin(2 + ( + u) sin(2 - ( - u) + — ^ sinr], (3.3) 
^ ^ 2 cos 1] 

sin(/7 — 2u] 

X2(u) = sin(2 + ( - u) sin(2 - ( - u) + fip, — sin{r] + 2m), 

2 cos rj 

/X •/« > \-/ri > \ ^sin(r7 — 2u) . 

Xsiu) = smU + C - m) smU + /^yW — sm r], 

^ ^ 2 cos rj 



where (,fi,fi, are arbitrary constants. If we choose p^{u) = p±{u) with 

^^^""^ " sin(2 + C ± w) sin(2 - C ± m) + ppsm'^'q ± 2u)/{2 cost]) ' ^^'^^ 
then K satisfies 

K{u)K{-u) = 1, K{0) = 1. (3.5) 

Class (B): sin2r7 ^ and sinAi] = 0. 
For this case we have three solutions: 

Xi{u) = sin(2 + ( + u) sin(| — — u), X2{u) = sin(2 + ( — u) sin(2 — ^ — m), 
X3{u) = sin(f + C-u) sin(§ -( + «), yi{u) = y2{u) = yi{u) = y2{u) = 0, (3.6) 
z{u) = z/ sin 2m, z{u) = 0, or z{u) = 0, z{u) = z/ sin 2m, 
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with z/ 7^ and ( (sin 2( ^ 0) being arbitrary constants, and 

yiiu) = y2{.u) = yi{u) = y2{u) = X2{u) = 0, 

xi{u) = 1, X3(m) = -1, z{u) = u, z{u) = -i, (3.7) 

where i/ 7^ is an arbitrary constant. 

If we take p^{u) = l/(sin(| + ( ± u) sin(| — C i 'u)) (four choices), then K-matrices 
( |3.6| ) satisfy ( |3.5| ). However, K- matrix ( |3.7|) does not satisfy (P-Sj). 



Class (C): sin2r7 = 0. 

For r/ = 0, TT, there are no restriction for K. For rj = ^vr, |7r, we have yi = y2 = fji 
y2 = and other components are arbitrary. 



4 XXZ Spin-1 Chain with General Boundary 

In this section we consider the K- matrix (|3.3|), K{u) = K{u; (, p, fi), and take p{u) = 
p±{u) and p^{u) = p±{u). K±{u) are defined as K^{u) = K{u; p-, p-), K^{u) = 
*ir(-M-r/;-C+,/i+,/i+). 

R- and K-matrices induce the hamiltonian of integrable spin-1 XX Z open chain with 
N sites, 

u=o trir+(0) ' ^ - ' 



N-l 



1 d 



n=l ^ "''^ 

where the two-site hamiltonian is given by 



Hn,n+l — ~;~Rn,n+l{u) 

du 



(4.2) 

u=0 



Here suffices of R and K± etc. indicate the sites on which the operators act. We remark 
that the terms in H which are not proportional to identity do not depend on the four 
choices (p±,p^). From (|;2D,(U),(pD,(|^, we obtain 

H = - — ^{sn- Sn+l — {Sn ■ 4+l)^ + ~ COS //) | S^'Sn+l ) 4 -^n+l + ^n^n+l} 

sm /?7 ^^-^ \ 

-(1 - cos 2r^) (44+1 - (sim^^r + (sir + (4+i)^)) 

(4.3) 



2 sin(| + C_) sin(2 — (_) + P-P- sin r/ tan 
X ((1 — P-P-) sinr/ s\s\ — sin2C_ s\ + \ ta.'n.ri{p't^s]^Si + p^^s'^s^) 

-y2sin(2 + C_)(/i„s+4 + f^-sl^i) - V2sin(2 - C_)(/i_4s+ + /i_s^4 
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~'~2 sin(| + C+) sin(| — C+) + sin tan 77 

X ((1 - sinr/ s^^rS^ - sin2C+ 3% + ^ tan?7(/i+s^s^ + /i+s]vSjv) 

-V2sin(f + C+)(/i+SivS?^ + fi+s%s^) - V2sin(f - C+)(/^+s^s^ + fi+s^s^ 
+ (constant) ■ id, 

where the spin-1 operator s^, s^{= ± is'^) is given by 



N 





( 1 





^ 






(0 


1 








/ 


o\ 























1 




= V2 


1 







lo 





-1 ) 






lo 





) 






[0 1 


0/ 



(4.4) 



In the special case, fi± = fi± = 0, ( [4. 3D reduces to the hamiltonian obtained in [TT| by 
using the diagonal K-matrix. 



5 Fusion Procedure 



The R-matrix ( p.2| ) of the 19-vertex model was obtained by directly solving the YBE in 



12|. It can also be constructed out of the R-matrix of the six- vertex model by using the 



fusion method 113 1. 

It was suggested that the fusion method can also be used to construct K-matrices of 
vertex models corresponding to higher spins out of the K-matrix of the six-vertex model 
In this section we apply this fusion method to the present case and compare the 
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result derived this way with that obtained above by directly solving the BYBE. 
The R-matrix associated with the spin-i representation of Uq{sl2) reads. 



RiU) 



( sin(?7 + u) 



V 





sinw 
sin?7 






sinr] 
sinw 










(5.1) 



sin(r7 + u) ) 



with pi{u) ^ being an arbitrary function. The corresponding reflection K-matrix has 



the form 



(1) K [ sin(^ + M) 1/ sin 2m 



z/ sin 2m sin(^ — u 

where ^, v are constant parameters and Pi{u) ^ is an arbitrary overall function. 



(5.2) 



According to |Tl[] the K-matrix for the 19-vertex model is given by the following fusion 
equation: 



K{u) = f{u)UP+ [k^'2)(u - f ) ® 1) R^'^'^\2u) (1 ® K^^\u + f )) P+U^\ 



(5.3) 
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where P+ is the projection onto the spin-1 component and U is a basis-changing matrix; 
exphcitly, 



P) 



/ 1 \ 

I I 
olio 
V 1 ; 



u 



/ 1 \ 

a a 

1 

V a -a J 



(5.4) 



It turns out that the K-matrix K{u) from (|5.3| ) agrees with one of the solutions in section 
3, eq.(^]3|), if one makes the following identification: 



sm{r] + 2u) pf (u - f )pf (u + f )pi (2m) ^ 



^ = C, ^= ,^ ^ = ,^ a = - — . 

y/ 2 COST] 2V 



(5.5) 



yJ2 COS ?7 ' 

It seems that the solutions in Class (B) can not be obtained from this fusion procedure. 
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A Derivation of the Reflection K-matrices 

In this Appendix we sketch how ( |3.3|) , ( p.6|) and ( p.7|) are obtained and show that there is 
no other nontrivial solution (Class (C) is trivial!). Our derivation is based on a careful 
case by case study. We sometimes omit explicit expressions when they may be lengthy. 
First we remark that one can check ( p.3| ), ( p.6| ) and (|3.7|) by direct substitution. Let us 
denote (i,j) component of (|3.1| ) as Eq[i,j]. Note that Eq[10 — i, 10 — j] can be obtained 



from Eq[2,j] by interchanging x\ ^ X3, y\ ^ ?/2, I/2 ^ Vi, z z. In the following we 
assume sin2?7 7^ (as mentioned above, the sm2ri = case is trivial). 

I. From Eq[2,8], we have 
sm{u-u')(yi{u)yi{u')-y2{u)y2{u')) = sm{ri -u-u')(yi{u')y2{u) -yi{u)y2{u')). (A.l) 



Dividing this by sin(M — u') sin^r] — u — u')y2{u)y2{u') and taking u' ^ u limit, we obtain 
the differential equation 

du 1 1/2 ) sin(?7 — 2u) 
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Solving this, we get 

yiju) _ sin(C- § + m) 

y^iu) sin(C + f-«)' ^ ■ ^ 

where ^ is a constant. We can check that this relation satisfies ( [A.l|) . 
Eq[2,9] gives rise to 

sin(r7 — u — u')y2{u)z{u') — sin(?7 — 2u')y2{u')z{u) = sm{u — u)yi{u)z{u'). (A. 3) 

By differentiating this linear equation with respect to u' and taking u' ^ u limit, we have 

1 dy2 1 dz\ 



yi = - cos 



{r] - 2u)y2 + sin(r7 - 2u)y2{- 



-y2 du z du' 

Solving this differential equation with the help of ( |A.2| ), we obtain 

y2{u) _ sin(C + f - 



Again we can check that this relation satisfies ( [A.3|) . Up to now, we have assumed y2 
and z are not zero. Taking into account the fact that they may be zero, we obtain the 
following result from ( [A.!] ) and ( |A.3D , 

yi{u) = /isin(C-f +M)fi'(M), y2{u) = /i sin(C+f z{u) = usm{T]-2u)g{u), (A.5) 

where (, fi, v are arbitrary constants and giu) ^ is an arbitrary function. 
Similarly, by Eq[8,2] and Eq[9,2], 

y^{u) = jx' sin(C - f + u)g{u), y2{u) = fi' sin(C + f - u)g{u), z{u) = v sin(?7 - 2u)g{u), 

(A.6) 

where (, fi', v' are arbitrary constants and g{u) ^ is an arbitrary function. 
From Eq[l,l] and Eq[9,9], we have 

sin(?7 + U + u')(yi{u')yi{u) - yi{u)yi{u')) = sin r](^z{u')z{u) - (A. 7) 

sin{'r] + u + u')(jj2{u')y2{u) - y2{u)y2{u)'j = sin rj (^z {u )z{u) - z{u)z{u')y (A. 8) 

Using (|A.5|) and ( [A.6|) , we arrive at the following three partial results: 



(A.9) 





' yi{u) 


= /Usin(C - ^ + u)g{u), yi{ 


u) = sin{( — ^ + u)g{u 




y2{u) 


= yUsin(C + f - u)g{u), y2{ 


u) = fl sin{( + ^ — u)g{u 






= u sin{ri — 2u)g{u), z{ 


u) = u sin^T] — 2u)g{u), 




' 1/1 (m) 


= sin(C - f +u)g{u), yi{u 


) = 0, 


(2) < 




= sin(C + f - u)g{u), y2{u 


) = 0, 




. An) 


= 0, z{u 


) = sin(?7 - 2m)^(m), 






= 0, yi{u) = 


= sin(C- ^ + u)g{u), 


(3) < 


y2{u) 


= 0, y2iu) = 


-- sin(C + f - u)g{u), 






= sin{ri — 2u)g{u), z{u) = 


-- 0, 



(A.IO) 



(A.ll) 
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where (, (, /i, fi, u, u are arbitrary constants and g{u) and g{u) are arbitrary functions 
{9:9^ 0, g/g ^ constant). 

II. First let us consider case (2). Eq[4,l] and Eq[2,l] can be written as the following 
form respectively 

Ai{u,u')g{u)g{u') + A2{u,u')g{u)g{u) = 0, 
A3{u,u')g{u)g{u') + A4{u,u')g{u')g{u) = 0, 

where Aj's are known functions of u and u' with their explicit expressions omitted. Since 
g{u)g{u') ^ 0, we have A{u,u') = A1A4 - A^A-^ = 0. A(0,77) = means sin2C = 0. 
With this however, A(w, —u) does not vanish, which implies the above two equations 
are contradictory each other. Therefore case (2) is not a solution. 

Similarly, from Eq[l,4] and Eq[l,2], we can show that case (3) is not a solution. 

III. In the following we consider case (1). Eq[l,4] is linear in x\{u)^X\{v!^^X2{u)^X2{u'^- 
There are two cases to be considered: // = and 7^ 0. For the former case, Eq[l,4] is 
satisfied if and only if \xv — 0. For the latter case, the coefficient of X2{v!) vanishes and 
x^iu) can be expressed by X\{u) and X\{u'^. Of course this expression must be independent 
on Vb . Similar results can be obtained for Eq[4,l]. Combining these results, we derive 

(i) ;U = 0, /i = 

(ii) = 0, /i ^ 1/ = 0) 

(iii) //^0,/i = (^i> = 0) (A.12) 
(iv-a) ^ 0, /i 7^ 0, 1/ = z> = 0) 

(iv-b) ;x^0,/i7^ 0,1/7^0 (^=^i>= (^)^, sin(7yT2C)7^o) 

(the condition sin(7^ + 2^) is obtained below). For (i) there is no restriction on xi, x^- 
For other cases, 

/ \ ( A2 \ . . , / , -/ii/ sin?7sin2C , , /* 

xi M = ^ + sm ^ „ ^g{u), A.13 

Vsmw cosM/ ^ [I sm(?7 — 2Q 

X2{u) = [— — )sin(f -C--»)^(^) + (^^f ■ ^ -sin2^)ff(^), 

where Ai, A2 arc arbitrary constants, and 5=1 for case (iv-b) and 5 = for other cases. 

Similarly from Eq[9,6] and Eq[6,9] we obtain ^3 and X2 (and the condition sin(?7+2C) 7^ 
in the case (iv-b) above): for (i) there is no restriction on 0:3, X2, and for other cases, 

/ \ / ^3 A4 \ . .„ ^ ... ./iz/ sinr/sin2C , , r , ■, a\ 

^3 = — ^ + Ml + C - u)g{u) - S——Y-—7^9{u), A.14 

\sm-u cos-u/ u smm + 2C 



sm{r] + 20' 

/ \ / -^3 -^4 \ . /„ A \ / \ rA'^ / sin^?7 . A / \ 

X2{u) = [- ]sm{^ + C-u)g{u) + 6—[————-sm2u)g{u), 

\smu cosuj ^ \i Vsm(?7-|-20 / 
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where A3, A4 are arbitrary constants, and S = 1 for case (iv-b) and 6 = for other cases. 
Combining (|A.13| ) and (|A.14|) (for example taking the hmit n ^ 0, ^), we get 
2A 

^li^) = ^^sin(f + C + M)sin(2 -( 
sm 2u 

.jiv sin?7 /^•„o/- sin(| — C — m) sin2?7 



H sm[7] — 2Q \ 2cosMsm(^ + Cj / 

2A 

X2{u) = ^— -sin(f + C-M)sin(f -C-u)^(u) (A.15) 
sm 

fiu sinr/ f . sin(f - C - sin2r/ N^ ■ o ( \ 

+V sin(,-2C) l""^+ 2cos.sin(f + C) j^H " 



2A 

^3(m) = ^^sm(2 -( + M)sm(f + C-M)fi'(M) 
sm2M ^ 



./iz/ sinr] ^ . sin(| + ^ — m) sin2r7 



+V sin(, + 2C) r - 2cos.sin(f-C) 

where 5 = 1 for case (iv-b) and 5 = for cases (ii),(iii),(iv-a). For case (i) there are no 
restriction for Xi, X2, 0:3. 

IV. Case (i). 

Let Xiiu) = Xi{u) / {sm{ri — 2u)g{u)). Eq[4,6] and Eq[6,4] lead to 

z/(Xi(m) - X3(iz) - (Xi(m') - X^{u'))) = 0, (A.16) 

V(x^{u) - X^{u) - (Xi(m') - X^{u'))) = 0, (A.17) 

and Eq[2,4] and Eq[6,8] to 

Xi(m)Xi(m') - X2{u)X2{u') + vv Xi{u)X2{u') - Xi(n')X2(u) 



sin(u + u') sm{u — u') 

X^{u)X^{u') - X2(n)X2(n') + vv X^{u)X2{u') - X;{u')X2{u) 



(A.18) 
(A.19) 



sin(u + u') sm{u — u') 

We have two choices (a) vi) = 0, (b) uu 7^ 0. 
IV. 1. Case (a). 

Through a similar calculation as for ( [A.l| ), it follows from (|A.18| ) and ( |A.19| ) that for 
X2 = 0, we have xi{u) = X2{u) = x^^u) = 0, and for X2 ^ 0, we get 

xi{u) = sin(ci — m) sin(c2 + m)/(m), 

X2{u) = sin(ci + u) sin(c2 + m)/(u), (A. 20) 

Xsi^u) = sin(ci + m) sin(c2 — «)/(«), 

where Ci, C2 are arbitrary constants and f{u) ^ is an arbitrary function. For the former 
case, we can check that they give the following two solutions of ( |3.1| ): 

z{u) = h{u), xi = X2 = x^i = yi = y2 = yi = m = z = 0, (A.21) 
z{u) = h{u), xi = X2 = X3 = yi = y2 = yi = y2 = z = 0, (A.22) 
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where h{u) ^ is an arbitrary function. For the latter case, Eq[3,5] imphes sin(ci + 
C2 + v) = 0- Since Xj's are invariant under C2 — ^ C2 + vr and f{u) —f{u), we can take 
Ci + C2 + rj = 0. Therefore we have 

xi{u) = sin(2 + ( + -u) sin(| - C - «)/(«), 

X2{u) = sin(2 + C-u) sin(f - C - u)f{u), (A.23) 
X3{u) = sin(| + C - m) sin(| - C + «)/(«), 

where ( is an arbitrary constant. The following three cases are possible: 

(al) u = i) = 0, 

(a2) z/ ^ 0, z> = 0, (A.24) 
(a3) u = 0, 9^0. 

• Case (al). We can check that (|3.1| ) is satisfied. This is a diagonal solution, 

yiiu) = y2{u) = yi{u) = y2iu) = z{u) = z{u) = 0, 
Xi{u) = sin(2 + ( + u) sin(2 — ( — u)f{u), 

X2iu) = sin(f + C - ^) sin(f - C - u)f{u), (A.25) 
X3{u) = sin(2 + C-u) sin(f - C + u)f{u), 

where ( is an arbitrary constant and f{u) ^ is an arbitrary function. 

• Case (a2). From ( |A. 16| ) we have Xi{u) —x^i^u) = a sin{ri — 2u)g{u), where a is a constant. 
From ( |A.23|) we get — sin 2^ sin 2m/(m) = crsin(?7 — 2u)g{u). There are two choices: 

(a2-l) sin2C 7^ (T ^ 0), 

(a2-2) sin2C = a = 0). (A.26) 

For the case (a2-l), Eq[l,3] implies sin4?7 = and then we can check that ( p.l|) is satisfied. 
Therefore we get a solution (note: sin 47] = 0), 

z{u) = z/ sin 2m h{u), yi{u) = y2{u) = yi{u) = y2{u) = z{u) = 0, 
xi{u) = sin(2 + ( + u) sin(2 - ( - u)h{u), 

X2{u) = sin(f + C-u) sin(f - C - u)h{u), (A.27) 
xs{u) = sin(2 + — u) sin(| — C + u)h{u), 

where z/ 7^ and ( (sin2C 7^ 0) are arbitrary constants and h{u) ^ is an arbitrary 
function. 

For the case (a2-2), we may take C = or ^tt because ( |A.23| ) is periodic in ( with 
period vr. For both the ( values, Eq[2,6] implies sm{ri — 2u)g{u) = (constant) x sin2'u/(-u). 
However Eq[l,3] implies this constant is zero, which contradicts g{u) ^ 0. 
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• Case (a3). From (|A.17| ) we have xi{u) — X3{u) = (Tsm{ri — 2u)g{u), where a is a 
constant. From ([A .231) we get — sin 2^ sin 2M/(n) = a sin{ri — 2u)g{u). Again two choices 
are possible: 

(a3-l) sin2C ^ a 7^ 0), 

(a3-2) sin2C = a = 0). (A.28) 

For the case (a2-l), Eq[3,l] imphes sin4?7 = and then we can check that is satisfied. 
Therefore we get a solution (note again: sin 4// = 0), 

z{u) = z>sin2M h{u), yi{u) = y2{u) = yi{u) = y2{u) = z{u) = 0, 
Xi{u) = sin(| + C + u) sin(2 - ( - u)h{u), 

X2{u) = sin(f + C-u) sin(f - C - u)h{u), (A.29) 
Xs{u) = sin(2 + ( — u) sin(| — C + u)h{u), 

where z/ 7^ and C, (sin2C 7^ 0) are arbitrary constants and h{u) ^ is an arbitrary 
function. For the case (a3-2), Eq[6,2] implies smij] — 2u)g{u) = (constant) x sin2u/(M). 
However Eq[3,l] implies this constant is zero, which contradicts g{u) ^ 0. 

IV. 2. Case (b). 

(|A.18| ) with u' = implies that Xi{u) + X2{u) is a constant, and (|A.18|) with u' = ^vr 
implies that Xi{u) — X2{u) is also a constant. Therefore Xi and X2 are both constants, 
Xi{u) = Ci, X2{u) = C2 with — C2 + z/z/ = 0. Similarly from ([A.19|) we have X^lu) = C3 
with cl — cl + uu = 0. We have two possibilities, C3 = Ci and C3 = — Ci. We eliminate z/ 
using = (c2 — Ci)/^. 

• Case C3 = ci. Eq[2,6] implies ci = and Eq[l,3] implies C2 = 0. But this contradicts 
z> ^ 0. 

• Case C3 = — Ci. Eq[3,7] implies ci = or C2 = 0. For ci = 0, Eq[l,3] implies C2 = 
but this contradicts z/ 7^ 0. For C2 = 0, Eq[l,3] implies ci = or sin4?7 = 0. The former 
contradicts z> 7^ 0. For the latter case we can check ( p.l|) is satisfied. Therefore we get a 
solution (for sin4?7 = 0), 

yiiu) = y2{u) = yiiu) = y2{u) = X2{u) = 0, 

xi{u) = h{u), X3{u) = —h{u), z{u) = phiu), z{u) = —^h{u)^ (A. 30) 
where z/ 7^ is an arbitrary constant and h{u) ^ is an arbitrary function. 

V. Case (ii), (in) 

First let us consider case (ii). Eq[6,2] implies fj? = AXucostj. Then we can check that 
( |3.1| ) is satisfied. Therefore we get a solution, 

yiiu) = y2{u) = z{u) = 0, 
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yi{u) = /isin(C - f +u)g{u), y2{u) = /ism(C + f - u)g{u), z{u) = usm{r] - 2u)g{u), 
xi{u) = ^\ sm(f + C + U) sin(f - C - u)g{u), (A.31) 

a;2(M) = 77^ sm(f + C - m) sin(2 - ( - u)g{u), 

2z/cos?7sm2M ^ 

X'i{u) = — — - sin(f + C -u) sin(f - C + u)g{u), 

Zu cos 7] sm 2u 

where v are arbitrary constant (/i 7^ 0, z/ 7^ 0) and g[u) ^ is an arbitrary function. 
Case (iii) is treated similarly: Eq[2,6] implies ^ = AXu cos rj and we get a solution, 

yi{u) = ij, sin(C - § + u)g{u), y2{u) = /i sin(C + f - u)g{u), z{u) = v sin(?7 - 1u)g{u), 
yi{u) = miu) = z{u) = 0, 

xi{u) = sin(f + C + U) sin(f - C - u)g{u), (A.32) 

2z/cos?7sm2M i i 

^^[u) = — - sin(f + C-u) sin(f - C - u)g{u), 

2z/cos?7sm2M ^ ^ 

3^3(m) = 7^ sin(2 + C - m) sin(2 -( + u)g{u), 

Zu cos 7] sm 2u 

where (, /i, u are arbitrary constant (/i 7^ 0, z/ 7^ 0) and (^('u) ^ is an arbitrary function. 

VI. Case (iv-a) 

Eq[l,7] implies /i = and Eq[7,l] implies /i = 0. But this contradicts the assumption. 
Therefore there are no solutions for this case. 

VII. Case (iv-b) 
Eq[2,4] implies 

^ ^ /x^ sin(f + C) sin(f - C) + f^J^^ sin rj sin 2r] 
4/iz/cos?7 sin(| + C) sin(| - C) 

Then we can check that is satisfied. Therefore we get a solution, 

yi{u) = IJ, sin(C - f + u)g{u), y2{u) = /i sin(C + f - M)fi'(M), ^(m) = sm{r] - 2u)g{u), 
yi{u) = /isin(C - f +u)g{u), y2{.u) = /isin(C + f - u)g{u), z{u) = z>sin(?7 - 2u)g{u), 

li^sm{^ + C + u)sm{^-C-u)+flu^sm2r]sm{r]-2u) r x oa\ 

xi{u) = ^ ^ — g{u), [AM) 

2/iz/ cos?7 sm 2u 

/i^ sin(| + ^ — m) sin(| — ^ — m) + 2/iz/^ cos?7 sin(?7 — 2-u) sin(77 + 2-u) 

X2[u) = — g{u), 

2/iz/ cos ?7sm 2u 

/i'^ sin(| + C — m) sin(| — C + n) + /2z/^ sin 2?7 sin(r7 — 2n) 

X3{u) = 7: ^ g{u), 

2jjLV COS?] sm 2u 

where v are arbitrary constants (/i 7^ 0, /i 7^ 0, z/ 7^ 0, z> = pP'u/ ji^) and g{u) ^ 

is an arbitrary function. 
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VIII. To summarize, we have obtained all the solutions of BYBE. ([OTD , ( 1091 ), (lOOl) 
correspond to three solutions in Class (B). Other solutions ( |A.21|) , ( |A.22|) , (|A.25|) , ( |A.31|) , 
( |A.32| ), ( [A.34| ) can be expressed in one formula ( p.3| ) (Class (A)). For example, ( |A.34| ) 



is obtained from (|3.3| ) by replacing ^ 2z//i ^ cos 77, fi — * 2/iz//i ^ cos 77 and dividing by 
2z//i~^ COST] sin 2-u. ( [A. 211 ) is obtained from ( ^.31 ) by setting /i = 0, dividing by /i^, and 



taking a limit — 00. But this solution can not satisfy ( p.5| ). 
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